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ABSTRACT 

Thermo chemical equilibrium and kinetic calculations are used to model sulfur 
and phosphorus chemistry in giant planets, brown dwarfs, and extrasolar giant 
planets (EGPs). The chemical behavior of individual S- and P-bearing gases and 
condensates is determined as a function of pressure, temperature, and metallicity. 
The results are independent of particular model atmospheres and, in principle, 
the equilibrium composition along the pressure-temperature profile of any ob- 
ject can be determined. Hydrogen sulfide (H 2 S) is the dominant S-bearing gas 
throughout substellar atmospheres and approximately represents the atmospheric 
sulfur inventory. Silicon sulfide (SiS) is a potential tracer of weather in substel- 
lar atmospheres. Disequilibrium abundances of phosphine (PH 3 ) approximately 
representative of the total atmospheric phosphorus inventory are expected to be 
mixed upward into the observable atmospheres of giant planets and T dwarfs. In 
hotter objects, several P-bearing gases (e.g., P 2 , PH 3 , PH 2 , PH, HCP) become 
increasingly important at high temperatures. 

Subject headings: astrochemistry — planets and satellites: individual (Jupiter) 
- stars: low-mass, brown dwarfs — stars: individual (Gliese 229B, HD 209458) 



1. Introduction 

The relatively cool, dense atmospheres of substellar objects - gas giant planets, brown 
dwarfs, and extrasolar giant planets (EGPs) - are ideal environments for the formation of 
molecules and condensates and the establishment of equilibrium chemistry. For this reason, 
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thermochemical models have been essential for interpreting and guiding spectroscopic obser- 
vations of the atmospheres of giant planets (e.g., see Lewis 1969a,b; Barshay & Lewis 1978; 
Fegley & Lewis 1979; Fegley & Prinn 1985, 1988; Fegley et al. 1991; Fegley & Lodders 1994) 
and brown dwarfs (e.g., Fegley & Lodders 1996). These models have also been useful for 
explaining how the optical and infrared spectra of substellar objects are influenced by CNO 
(Lodders & Fegley 2002), alkali (Lodders 1999a), Ti and V (Lodders 2002), and condensation 
(e.g., Lodders & Fegley 2006) chemistry. 

After H, C, N, and O, sulfur and phosphorus are the next most abundant chemically 
reactive volatile elements in a solar system composition gas. Fegley & Lodders (1994; here- 
after FL94) modeled the chemistry of all naturally occurring elements, including S and P, 
along the atmospheric profiles of Jupiter and Saturn. Recently, Lodders & Fegley (2002; 
hereafter LF02) modeled CNO chemistry in substellar objects as a function of pressure, 
temperature, and metallicity. Here we continue and extend these studies by employing ther- 
mochemical equilibrium and kinetic calculations to model sulfur and phosphorus chemistry 
in the atmospheres of substellar objects. 

Our approach is similar to that of LF02. We compute the chemistry of S- and P-bearing 
gases and condensates as a function of total pressure, temperature, and metallicity, and our 
abundance results are independent of any particular atmospheric profile. In principle, the 
pressure-temperature profile for any substellar object may be superimposed on our abun- 
dance contour diagrams to find the equilibrium composition along the profile. 

The paper is organized as follows. In section (§2) we describe our computational meth- 
ods. In §3 we describe sulfur chemistry in a solar system composition gas: we identify 
major gases and condensates, and discuss the chemistry of individual S-bearing species as a 
function of total pressure, temperature, and metallicity. We then apply our results by exam- 
ining the sulfur species along the atmospheric profiles of representative substellar objects. 
A similar treatment of phosphorus chemistry follows in §4. We discuss the relative roles of 
thermochemistry and photochemistry in §5, and conclude with a brief summary (§6). 

2. Computational Methods 

Thermochemical equilibrium calculations were performed using a Gibbs free energy 
minimization code, previously used for modeling the atmospheric chemistry of Saturn (Viss- 
cher & Fegley 2005). Where relevant, we considered the effects of vertical mixing on the 
abundances of gases (e.g., PH 3 ) which serve as chemical probes of the deep atmospheres of 
Jupiter and Saturn. This was done using a chemical dynamical model described in Fegley 
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& Prinn (1985). We generally focus on temperatures of 1000 K and higher, where thermo- 
chemistry is expected to be much more important than photochemistry in the atmospheres 
of extrasolar giant planets close to their primary stars. Further discussion of thermochemical 
vs. photochemical processes is given in §5. 

Thermodynamic data for the equilibrium calculations were taken from the compilations 
of Gurvich et al. (1989-1994), Robie & Hemingway (1995), the fourth edition of the JANAF 
Tables (Chase 1999), and the thermodynamic database maintained in the Planetary Chem- 
istry Laboratory (FL94; LF02). This database includes additional thermodynamic data 
from the literature for compounds absent from the other compilations, as well as several 
important data revisions for the sulfur and phosphorus species SH, S2O, NS, PS, PH, PH 3 , 
PN, and Mg 3 P 2 8 (s), which are incorrect in the JANAF tables (Lodders 1999b, 2004a). 
Thermodynamic data for P 4 6 was taken from the JANAF tables for reasons cited in FL94. 

All calculations were done with solar system (i.e., protosolar) elemental abundances from 
Lodders (2003). The effect of metallicity on sulfur and phosphorus chemistry was examined 
by running computations at [M/H] = -0.5 dex (subsolar), [M/H] = dex (solar), and [M/H] 
= +0.5 dex (enhanced) metallicities, where M is any element of interest (e.g., S, P, C, O) 1 . 
The metallicity factor, m, is defined as logm = [M/H]. Stellar abundance determinations 
show that element abundance ratios [M/H] for elements of interest here vary similarly with 
[Fe/H] between -0.5 < [Fe/H] < +0.5 dex (Edvardsson et al. 1993; Gustafsson et al. 1999; 
Chen et al. 2000; Smith et al. 2001; Ecuvillon et al. 2004; Huang et al. 2005), consistent 
with what one may expect from galactic chemical evolution models (e.g., Timmes et al. 
1995). The slopes for fits of [M/H] vs. [Fe/H] are about 0.61, 0.65, and 0.6 for M = C, 
S, and O, respectively, which justifies our approach of "uniform" metallicity variations for 
these elements (e.g., [S/H] « [P/H] [C/H] « [O/H] = [M/H]) over the metallicity range 
considered here. Because phosphorus has weak lines, galactic chemical evolution models are 
used to estimate P abundances as a function of metallicity (e.g., Timmes et al. 1995; Samland 
1998). We thus assume a similar enrichment in phosphorus as for other heavy elements. 
However, the effects of metallicity on chemical equilibrium abundances are explicitly treated 
below and can be used to study effects from non-uniform elemental abundance enrichments. 

For each S- and P-bearing species, we discuss how the most plausible representative net 
thermochemical formation reaction depends on temperature, total pressure, and metallicity. 
For example, consider the formation of thioxophosphino (PS) via the net thermochemical 
reaction 

PH 3 + H 2 S = PS + 2.5H 2 



1 We use the conventional notation [M/H] = log(M/H) - log(M/H) solar . 
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This endothermic reaction proceeds toward the right with increasing temperature and the 
PS abundance increases. According to LeChatelier's principle, this reaction also proceeds 
toward the right with decreasing pressure because there are 2 gas molecules on the left and 
3.5 gas molecules on the right. Writing out the equilibrium constant (K p ) expression for PS 
in terms of the total pressure (Pt) and mole fractions 2 

Xps = (Xp^X^s/X^KpP^ 1 - 5 , 

shows that the mole fraction of PS is proportional to P^ 1 ' 5 . The metallicity dependence 
of Xps enters only through the PH 3 and H 2 S abundances, whereas the H 2 abundance is 
metallicity-independent by definition. Phosphine and H 2 S are typically the most abundant 
P- and S-bearing gases, respectively. Hence Xpn 3 and Xh 2 s must each be expanded with a 
metallicity factor m, which results in an overall m 2 metallicity dependence for Aps, assuming 
[S/H] = [P/H] (see above). This basic approach is used to describe the equilibrium chemical 
behavior of each S- and P-bearing species throughout the paper. 



3. Sulfur Chemistry 

3.1. Sulfur Gas Chemistry 

The sulfur equilibrium gas chemistry as a function of temperature and total pressure in 
a protosolar composition gas is illustrated in Figure 1. The most abundant gases are H 2 S, 
SH, or monatomic S, and the P-T regions where each gas is dominant are bounded by solid 
lines. Also shown are condensation curves (dotted lines) for S-bearing condensates (see §3.2) 
and the H 2 =H and CH 4 =CO equal abundance curves (dash-dot lines). Model atmosphere 
profiles for Jupiter (T eff = 124 K, \ogg = 3.4), the T dwarf Gliese 229B (T eff = 960 K, 
\ogg = 5.0; Marley et al. 1996), the close-orbiting EGP (or "Pegasi" planet) HD209458b 
(T cff = 1350 K, \ogg = 3.0; Iro et al. 2005), an L dwarf (T cfr = 1800 K, \ogg = 5.0; Burrows 
et al. 2006), and an M dwarf (T cff = 2600 K, dust-free, \ogg = 5.0; Tsuji et al. 1996) are 
indicated by dashed lines. The overall chemistry in the atmospheres of Jupiter and Gliese 
229B differs slightly from the chemistry of a solar metallicity gas because Jupiter has a heavy 
element enrichment comparable to [M/H] +0.5 dex (Lodders 1999a; LF02), whereas Gliese 
229B likely has a subsolar metallicity of [M/H] « -0.3 dex (Saumon et al. 2000); their P-T 
profiles are shown here for reference. However, as described below, the P-T boundaries for 
the sulfur gases in Figure 1 are metallicity independent. This is because the P-T boundaries 



2 The mole fraction abundance for a gas i is defined as X t = total moTs'oVall gases " 
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are defined by equal abundance ratios of the neighboring gases (e.g., H 2 S/SH, SH/S, H 2 S/S), 
so that the metallicity dependence cancels out. 

Hydrogen sulfide (H 2 S) is the dominant S-bearing gas in substellar atmospheres under 
most P-T conditions considered here. The absolute H 2 S abundance is of course metallicity 
dependent (see §3.1.1). At high temperatures and low pressures, H 2 S is replaced by SH via 
the net thermochemical reaction 

H 2 S = SH + 0.5H 2 , (1) 

and SH becomes the dominant sulfur gas. The equilibrium constant expression for reaction 
(1) is 

K x = (XshX0 2 5 /X H2 s)P£- 5 . (2) 

Rearranging and substituting for the temperature dependence of K\ (logfTi = 3.37 — 8785/T 
from 800 to 2500 K), the H 2 S/SH ratio is given by 

log(X H2 s/Xs H ) = -3.37 + 8785/T + 0.5 logP T + 0.5 logX H2 . (3) 

The H 2 S/SH ratio is independent of metallicity because the m dependence of each S-bearing 
gas in reaction (1) cancels out and the H 2 abundance is essentially constant (X H2 ~ 0.84) 
over small metallicity variations. The solid line separating the H 2 S and SH fields in Figure 1 
shows where these two gases have equal abundances. With Xu 2S /X S u = 1 and Xu 2 ~ 0.84, 
equation (3) can be rewritten to give the H 2 S=SH equal abundance boundary: 

log P T = 6.82 - 17570 /T, (4) 

which is independent of metallicity. The abundance of each sulfur gas does not drop to zero 
as this line is crossed; H 2 S is still present within the SH field and vice versa (see §3.1.2). 

With increasing temperatures, SH dissociates and monatomic S becomes the dominant 
sulfur-bearing gas via 

SH = S + 0.5H 2 . (5) 

Using logics = 2.36 - 7261/T from 800 to 2500 K and rearranging, the SH/S abundance 
ratio is given by 

log(X SH /X s ) = -2.36 + 7261/T + 0.5 log P T + 0.5 log X H2 . (6) 

This ratio is also independent of metallicity because the m dependence of each S-bearing gas 
in reaction (5) cancels out. With Xsh/Xs = 1 and Xh 2 ~ 0.84, equation (6) is rewritten to 
give the position of the metallicity-independent SH=S equal abundance line 



logP T = 4.80 - 14522/T. 



(7) 
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At the lowest total pressures considered here, monatomic S is the major sulfur gas (see 
Figure 1). At temperatures below ~ 1509 K, H 2 S directly converts to monatomic S via the 
net thermochemical reaction 

H 2 S = S + H 2 . (8) 
The H 2 S/S abundance ratio follows from adding equations (3) and (6) as 

log(X H2S /X s ) = -5.73 + 16046/T + \ogP T + logX H2 , (9) 

and is independent of metallicity. Again, with Xn 2 s/Xs = 1 and Xh 2 ~ 0.84, the metallicity- 
independent position of the H 2 S=S boundary is approximated by 

log P T = 5.81 - 16046/T. (10) 

The H 2 S=SH, SH=S, and H 2 S=S boundaries intersect at the H 2 S-SH-S "triple point" at 
T ~ 1509 K and Pt ~ 10~ 4 ' 82 bar, indicated by the triangle in Figure 1 and represented by 
the net thermochemical equilibrium 

H 2 S + S = 2SH. (11) 

The P-T position of the triple point is independent of metallicity because the intersecting 
boundaries in equations (4), (7), and (10) are each independent of metallicity. At the triple 
point, all three gases have equal abundances A(H 2 S) = A(SH) = A(S) « |ES (where SS is 
the total elemental sulfur abundance) of one third of the total sulfur abundance. 

As noted above, the H 2 abundance is essentially constant (X H2 ~ 0.84) over the small 
metallicity range considered here. However, at high T and low P T) H 2 thermally dissociates 
to H via the net thermochemical reaction 

H 2 = 2H. (12) 

The position of the H 2 =H equal abundance boundary is indicated by the dash-dot line in 
Figure 1. The abundance of H 2 as a function of Pt and T in a solar system composition gas 
can be determined by solving the expression 

1.9845 + 10 M - v/10"(3.9690 + 10") 
XH2 * 2^670 ' (13) 

where the variable u is given by 

u = -23672/T - log P T + 6.2645. (14) 

At pressures and temperatures relevant for substellar atmospheres (Figure 1), 10" approaches 
zero and X H2 ~ 0.84. 
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3.1.1. Hydrogen Sulfide, H 2 S gas 

As shown in Figure 1, H2S is expected to be the most abundant S-bearing gas throughout 
the atmospheres of substellar objects (e.g. FL94; Fegley & Lodders 1996). Aside from minor 
amounts of sulfur (< 9%) removed by metal sulfide clouds (see §3.2), H 2 S approximately 
represents the atmospheric sulfur inventory until it is removed by very low temperature cloud 
formation and/or photochemical destruction. The H 2 S abundance also slightly decreases at 
atmospheric levels below the Mg-silicate clouds where SiS gas is relatively abundant (see 
§3.1.4 and §3.3). Within the H 2 S field and at the high temperatures below the metal sulfide 
clouds, the H 2 S abundance is given by 

logX H2 s~-4.52 + [S/H], (15) 

whereas above the metal sulfide clouds, the H 2 S abundance is 

logX H2S ~-4.56 + [S/H]. (16) 

Figure 2a gives mole fraction contours (on a logarithmic scale) of H 2 S as a function of 
temperature and total pressure in a solar metallicity gas. The H 2 S abundance at higher or 
lower metallicities can be found by substituting for [S/H] in equations (15) and (16). As 
discussed above, H 2 S dissociates to SH and S at high temperatures. Inside the SH and S 
fields, the H 2 S abundance remains proportional to m and decreases when moving toward 
higher temperatures and lower total pressures. 

3.1.2. Mercapto, SH gas 

Mole fraction contours of SH as a function of T and Pt are shown in Figure 2b. Within 
the H 2 S field (ATh 2 s ~ A^g and Xh 2 ~ 0.84), mercapto is governed by reaction (1). The SH 
abundance is given by substituting equation (15) into equation (3) and rearranging 

logXgH « -1.11 - 8785/T- 0.5 log P T + [S/H], (17) 

proportional to P T 5 and m. At high temperatures, SH is replaced by S via reaction (5). 
Within the S field, the SH abundance is proportional to m and decreases when moving away 
from the SH-S boundary to higher T and lower Pp. 

3.1.3. Monatomic Sulfur, S gas 

Mole fraction contours for monatomic S are shown in Figure 2c. Inside the H 2 S field, 
monatomic sulfur gas is governed by reaction (8). The S abundance is given by rearranging 
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equation (9), using logX H2 s from equation (15) and logX H2 ~ —0.08 to give 

logXg « 1.29 - 16046/T - logP T + [S/H], (18) 

proportional to P^ 1 and m. As illustrated in Figure 1, monatomic S is the dominant S- 
bearing gas (X$ & X^s) at high T and low Pt in a solar composition gas. Thermal ionization 
of S only becomes important at temperatures higher than those shown in Figure 1, e.g., 
S+/S ~ 1 at 6974 K (10~ 2 bars), 5340 K (10~ 4 bars), and 4326 K (10~ 6 bars). 



3.1.4- Silicon Sulfide, SiS gas 

The chemistry of SiS in a solar metallicity gas is shown in Figure 2d. Also shown are 
dashed lines for the condensation curves of forsterite (Mg 2 Si04) and enstatite (MgSiOs). 
Condensation of Mg-silicates efficiently removes Si gases from the atmospheres of gas giant 
planets and cool brown dwarfs (Fegley & Prinn 1988; FL94). The SiS abundance is controlled 
by the net thermochemical reaction 

SiO + H 2 S = SiS + H 2 0. (19) 

The SiS abundance at high temperatures below the Mg-silicate clouds is given by 

logX SiS w -1.07 + 666/T + [M/H] + logX H2 s, (20) 

assuming [O/H] = [S/H] (see §2). Below the clouds, X^is is relatively constant (~ 5 ppm at 
[M/H] = 0) until H 2 S is replaced by monatomic S at high T and low Pt- The SiS abundance 
also decreases at Pt > 100 bar as CO is replaced by CH 4 (see Figure 1). This conversion 
increases the H 2 abundance (see LF02), and drives reaction (19) toward the left to yield 
less SiS. 

At lower temperatures, above the clouds, the SiO abundance is governed by the reaction 

Mg 2 Si0 4 (s) + SiO + H 2 = 2MgSi0 3 (s) + H 2 . (21) 

Combined with reaction (19), the abundance of SiS is given by 

logXsis ~ 5.38 - 28151/T - logP T - 21ogX H2 o + [S/H]. (22) 

Although reaction (19) does not indicate a dependence on P T and suggests a linear depen- 
dence on the H 2 abundance, the SiS abundance in equation (22) is proportional to P T X 
and Ay 2 2 Q . These dependencies enter through the SiO abundance governed by reaction (21). 
The H 2 abundance above the Mg-silicate clouds is approximated by 



logX H2 o « -3.58 + 0.46/(1 + 10") + [O/H], 



(23) 
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where 

u = — 11704/T — 2 log Px + 9.78. (24) 

Substituting logX H2 o from equation (23) into equation (22) shows that X SiS is proportional 
to m -1 , assuming [O/H] = [S/H] (see §2). Curvature in the SiS abundance contours results 
from the dependence of the H2O abundance on the CH4/CO equilibrium (see LF02). For 
reference, X H2 o ~ 10~ 312 in CH 4 -dominated objects and X H2 o ~ 10~ 3 58 in CO-dominated 
objects with solar metallicity in atmospheric regions above the Mg-silicate clouds. 



3.2. Sulfur Condensation Chemistry 

With decreasing temperatures, sulfur condenses into MnS, Na 2 S, ZnS, and NH 4 SH cloud 
layers. Condensation affects the spectra of substellar objects by removing gases (and thus 
potential opacity sources) from the observable atmosphere, and by introducing cloud particles 
(Marley et al. 1996; Lodders 1999a; Burrows et al. 2000a; Lodders & Fegley 2006). 

The formation of a sulfur-bearing cloud requires that its constituent elements were not 
removed by condensation into more refractory cloud layers at deeper levels. For example, 
one might naively expect FeS cloud formation because FeS is a well-known condensate in 
the equilibrium condensation sequence of the solar nebula gas. However, condensation of 
FeS is preceded by Fe metal condensation and in the high-gravity environment of substellar 
objects, iron settles into a deep cloud layer at high temperatures (e.g. Prinn & Olaguer 1981; 
FL94; Lodders 1999a; LF02; Lodders & Fegley 2006). Hence no Fe metal is left at lower 
temperatures when FeS could condense. Note that if FeS condensation did occur, all H 2 S 
would be removed because Fe/S 1.9. However, the measured H 2 S abundance of ~ 2.4 
times the protosolar H 2 S/H 2 ratio in Jupiter's atmosphere (Niemann et al. 1998; Wong et 
al. 2004) requires Fe metal cloud formation at deep atmospheric levels. 

Sulfur thus condenses as other sulfides with metals (Mn, Na, Zn) that are not already 
removed into high-temperature condensates. The abundances of Mn, Zn, and Na are all 
lower than that of S and condensation of these sulfides only depletes the atmosphere of 
~ 9% of its entire sulfur inventory. Therefore, H 2 S gas remains in T dwarfs, L dwarfs, and 
Pegasi planets, and is only depleted in objects with low enough temperatures for NH4SH 
condensation. Table 1 lists the calculated condensation temperatures of MnS, Na 2 S, ZnS, 
and NH 4 SH for HD209458b and a model L dwarf with [M/H] « 0, for Jupiter ([M/H] « 
+0.5), and for Gliese 229B ([M/H] -0.3), assuming uniform enrichments and/or depletions 
of heavy element abundances. 
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3.2.1. Alabandite, MnS, condensation 

The most refractory S-bearing condensate in substellar atmospheres is MnS (Figure 1). 
Alabandite clouds form via the net thermochemical reaction 

Mn + H 2 S = MnS(s) + H 2 , (25) 

where the condensation temperature is approximated by 

10 4 /T cond (MnS) « 7.45 - 0.42(logP T + [Mn/H] + [S/H]) 

« 7.45 -0.42 log P T - 0.84[M/H], (26) 

assuming [Mn/H] = [S/H] (see §2). Because MnS condenses at sufficiently high temperatures 
(e.g. ~ 1340 K at 1 bar), MnS clouds are expected in most substellar atmospheres. The Mn 
abundance is 2% that of sulfur and thus MnS cloud formation removes all Mn and 2% of all 
sulfur from the atmosphere. 



3.2.2. Sodium Sulfide, NCL2S, condensation 

The Na 2 S cloud layer forms via the net thermochemical reaction 

2Na + H 2 S = Na 2 S(s) +H 2 , (27) 

The condensation temperature of Na 2 S as a function of Pt and metallicity is approximated 
by 

10 4 /T cond (Na 2 S) « 10.05 - 0.72(logP T + [Na/H] + 0.5[S/H]) 

« 10.05 - 0.72 log P T - 1.08[M/H], (28) 

assuming [Na/H] = [S/H] (see §2), and is shown in Figure 1 for a solar-metallicity gas. Below 
the Na 2 S cloud, monatomic Na is the most abundant Na-bearing gas, followed closely by 
NaCl. With decreasing temperatures, NaCl becomes increasingly important, but its abun- 
dance never exceeds that of monatomic Na below the cloud (Lodders 1999a). Sodium sulfide 
condensation effectively removes all sodium and ~ 6.5% of all sulfur from the atmosphere 
because the protosolar Na abundance is 13% of the protosolar S abundance. 



3.2.3. Sphalerite, ZnS, condensation 

Sphalerite clouds form by the net thermochemical reaction 

Zn + H 2 S = ZnS(s) + H 2 . 



(29) 
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The condensation temperature of ZnS as a function of Pt and metallicity is given by 

10 4 /^cond(ZnS) « 12.52 - 0.63(logP T + [Zn/H] + [S/H]) 

« 12.52 - 0.63 log P T - 1.26[M/H], (30) 

assuming [Zn/H] = [S/H] (see §2). Sphalerite cloud formation effectively removes all Zn and 
~ 0.3% of H 2 S from the atmosphere. 



3.2.4- Ammonium Hydrogen Sulfide, NH^SH, condensation 

Sulfur condenses as NH 4 SH via the net thermochemical reaction 

NH 3 + H 2 S = NH 4 SH(s), (31) 

in the cool upper atmospheres of giant planets (Figure 1), where H2S and NH3 are the major 
S- and N-bearing gases, respectively. The NH 4 SH condensation temperature as a function 
of Pt and metallicity is approximated by 

10 4 /^cond(NH 4 SH) « 48.91 - 4.15(logP T + 0.5 [N/H] + 0.5 [S/H]) 

« 48.91 -4.15 log P T -4.15[M/H], (32) 

assuming [N/H] = [S/H] (see §2). This expression differs slightly from LF02 due to the 
revised solar elemental abundances of Lodders (2003). The formation of a NH 4 SH cloud is 
expected to efficiently remove all sulfur from the atmosphere because NH3 is ~ 4.4 times 
more abundant than H 2 S in a protosolar composition gas. 



3.2.5. Hydrogen Sulfide, H 2 S, condensation 

Hydrogen sulfide cloud formation is only possible in cool objects (such as Uranus and 
Neptune) if NH 3 is absent or less abundant than H 2 S, since ammonia will consume atmo- 
spheric sulfur to form NH 4 SH via reaction (31) (e.g., see Fegley et al. 1991). At sufficiently 
low temperatures, H 2 S condensation may occur via the net reaction 

H 2 S = H 2 S (s). (33) 

The H 2 S condensation temperature is approximated by 

10 4 /^cond(H 2 S) « 86.49 - 8.54 log P T - 8.54[S/H]. (34) 

We found that no other S-bearing condensates are stable (OCS, S0 2 , elemental S) over this 
P- T range in a protosolar composition gas. 
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3.3. Sulfur Chemistry in Substellar Objects 

Figures 3a-3d summarize the equilibrium sulfur gas chemistry computed along the P-T 
profiles of four representative substellar objects: Jupiter (T e g = 124 K), the T dwarf Gliese 
229B (T eS = 960 K), the Pegasi planet HD209458b (T cff = 1350 K), and an L dwarf (T eff = 
1800 K) (see §3.1). The profiles for Jupiter and Gliese 229B were computed at [M/H] « +0.5 
and [M/H] —0.3, respectively, assuming uniform enrichments and depletions in elemental 
abundances. The upper atmosphere of HD209458b is heated by a large incoming stellar 
flux, which is responsible for the shape of this planet's P-T profile (Figure 1) and hence the 
curvature in all gas abundances near 1800 K . Also shown are gases which are strongly affected 
by S-bearing cloud formation (e.g., Na, NaCl, Mn, Zn). The condensation temperatures of 
Mg2Si04 (Fo), MgSiOs (En), MnS, Na2S, and ZnS are indicated by arrows (see Table 1). 
The Jupiter profile intersects the forsterite condensation curve at high temperatures outside 
the range of Figure 3a. The Gliese 229B profile intersects the ZnS condensation curve at 780 
K, outside the range of Figure 3b. In contrast, ZnS does not condense in the atmospheres 
of HD209458b or the L dwarf. 

Hydrogen sulfide is the most abundant S-bearing gas throughout the atmospheres of 
substellar objects (X H2 s ~ X-^s)- The H 2 S abundance slightly increases when silicon in SiS 
is removed by Mg-silicate cloud formation (e.g., see Figures 3b-3d), and slightly decreases 
when sulfur is removed by metal sulfide cloud formation (e.g., see Figure 3a). 

The Galileo entry probe measured H 2 S down to the ~ 20 bar level on Jupiter (Niemann 
ct al. 1998; Wong et al. 2004), where the observed H 2 S abundance of 2.4 times the solar 
H 2 S/H 2 ratio should represent the total sulfur inventory in Jupiter's atmosphere (e.g., Bar- 
shay & Lewis 1978; FL94; Lodders 2004b). Ground-based observations have failed to detect 
H 2 S on Jupiter and Saturn because of its removal by condensation, its short photochemical 
lifetime, and the lack of major H 2 S features in the 5 /im window. Observations of Jupiter at 
2.7 /im by KAO (Larson et al. 1984) and at 8.5 /xm by Voyager (Bezard et al. 1983) provide 
only upper limits because these wavelengths probe atmospheric levels where H 2 S is depleted 
by NH 4 SH condensation. The HST detection of H 2 S in the Shoemaker-Levy 9 impact region 
suggested excavation of Jovian atmosphere from regions below the NH 4 SH cloud where H 2 S 
is more abundant (Noll et al. 1995; Zahnle et al. 1995). 

Based on tropospheric NH 3 and H 2 S abundances, NH 4 SH cloud formation is expected 
to occur near the 220 K, 2.4 bar level on Jupiter (cf. FL94). This should efficiently remove 
H 2 S from the gas because NH 3 /H 2 S ~ 7.5 ± 3.4 in Jupiter's atmosphere (Wong et al. 2004). 
The inferred NH 4 SH cloud layer detected by the Galileo entry probe nephelometer was thin 
and tenuous, presumably due to the probe's entry in a relatively clear 5 /im hot spot (West 
et al. 2004). 
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In contrast, NH 4 SH cloud formation is not expected in the warmer atmospheres of brown 
dwarfs or Pegasi planets because their P-T profiles do not intersect the NH 4 SH condensation 
curve (see Figure 1), and H 2 S gas remains. Measurements of H 2 S in these objects should 
therefore approximately represent their bulk atmospheric sulfur inventory (e.g. Fegley & 
Lodders 1996). Saumon et al. (2000) investigated the observability of H2S and found that a 
H 2 S feature at 2.1 //m is potentially detectable in the infrared spectrum of Gliese 229B. In 
the uppermost atmospheres of brown dwarfs and Pegasi planets, H 2 S is plausibly destroyed 
by photochemistry (see §5), as is expected for Jupiter (e.g., Prinn & Owen 1976). 

Over the P-T range considered here, the next most abundant S-bearing gases after 
H 2 S are SH (at low temperatures) or SiS (at high temperatures). The mercapto radical is 
generally the second or third most abundant sulfur gas throughout substellar atmospheres. 
Yamamura et al. (2000) identified SH at ~ 3.67 fim in the published infrared spectrum 
of R Andromedae, a S-type star. Their inferred SH/H ratio of ~ 10~ 7 is consistent with 
thermochemical equilibrium (Tsuji 1973, and Figure 2b). Berdyugina & Livingston (2002) 
identified SH at ~ 3280 A in the Sun's photosphere. They do not give column densities, 
but computed synthetic spectra at 5250 and 4750 K. We calculated SH mole fractions of 
~ 9 x 10~ 9 (5250 K) and ~ 5 x 10~ 9 (4750 K) using the solar model atmosphere of Edvardsson 
et al. (1993). Using the effective temperatures for reference, we also compute SH column 
densities of ~ 7 x 10 12 (960 K), ~ 2 x 10 17 (1350 K), and - 2 x 10 17 cm" 2 (1800 K) on Gliese 
229B, HD209458b, and the L dwarf, respectively, versus the SH column density of ~ 4 x 10 20 
cm -2 in R Andromedae (Yamamura et al. 2000). The equilibrium SH abundances expected 
for brown dwarfs and Pegasi planets may be too low for reliable quantitative analysis. 

The relative importance of SiS increases with increasing temperature, and peak SiS 
abundances of ~ 6 ppm at [M/H] = and ~ 60 ppm at [M/H] = +0.5 are achieved at 
P- T conditions close to the condensation temperatures of the Mg-silicate clouds. With their 
condensation, silicon is removed from the atmosphere and the abundance of SiS and other 
Si-bearing gases drop. This makes SiS a potential tracer of weather in Pegasi planets and 
L dwarfs, analogous to FeH in T dwarfs (Burgasser et al. 2002), and gaps in the Mg-silicate 
cloud layers may expose regions where SiS is relatively abundant. In cooler objects, the SiS 
fundamental absorption band near 13.3 fim may be difficult to distinguish from ammonia 
features as NH 3 replaces N 2 to become the dominant N-bearing gas; the first SiS overtone 
at 6.6 /im has been observed in the carbon giant WZ Cas (Aoki et al. 1998), though this 
feature will likely be overwhelmed by H 2 absorption in substellar atmospheres. 

The increasing importance of NaCl instead of Na and the removal of Na by Na 2 S cloud 
formation is a plausible cause of the observed weakening of Na atomic lines throughout the 
L dwarf spectral sequence and their disappearance in early T dwarfs (e.g. Kirkpatrick et 
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al. 1999; McLean et al. 2003); this is further enhanced by pressure broadening of the Na I 
doublet (e.g. Tsuji et al. 1996; Burrows et al. 2000b; Liebert et al. 2000). Sodium sulfide is 
expected to condense at the ~ 1370 K level in Jupiter's atmosphere and the ~ 1000 K level 
in Gliese 229B (see Table 1). The L dwarf and the Pegasi planet HD209458b are too warm 
for Na2S condensation (assuming [M/H] ~ 0; see Figure 1) and Na remains in the gas. This 
is consistent with the detection of Na in the atmosphere of HD209458b (Charbonneau et 
al. 2002), though the chemistry results of Lodders (1999a) suggest that limited Na 2 S cloud 
formation may occur on this planet if night-side temperatures are low enough (e.g., Fortney et 
al. 2003, 2005; Iro et al. 2005; Barman et al. 2005), or if its metallicity is sufficiently enhanced. 
However, while Na chemistry is strongly affected by Na 2 S condensation, S chemistry remains 
almost unchanged because Na 2 S cloud formation only removes 6% of all sulfur (see §3.2.2). 

The formation of the MnS cloud layer effectively removes Mn from the atmospheres of 
substellar objects. Similarly, ZnS condensation removes Zn from the atmosphere of Jupiter 
([M/H] « +0.5) and Gliese 229B ([M/H] w -0.3) above the ~ 980 K and - 780 K levels, 
respectively. In contrast, the atmospheres of L dwarfs and Pegasi planets are too warm for 
ZnS condensation and Zn will remain in the gas. 



4. Phosphorus Chemistry 

4.1. Phosphorus Gas Chemistry 

The phosphorus equilibrium gas chemistry as a function of temperature and total pres- 
sure in a protosolar composition gas is illustrated in Figure 4. The P-T regions indicating 
the most abundant phosphorus bearing gas are bounded by solid lines. Also shown is the 
condensation curve for NH 4 H 2 P0 4 (dotted line), model atmosphere profiles for representa- 
tive substellar objects (dashed lines; see §3.1), and the H 2 =H and CH 4 =CO equal abundance 
curves (dash-dot lines). Phosphorus chemistry is considerably more complex than that for 
C, N, O, and S, and several P-bearing species become abundant at elevated temperatures 
(Barshay & Lewis 1978; Fegley & Lewis 1980; FL94). 

Under equilibrium conditions, P40 6 is the dominant phosphorus gas at low tempera- 
tures, and is replaced by PH 3 or P 2 as temperatures increase. The conversion of P 4 6 to 
PH 3 occurs by net thermochemical reaction: 

P 4 6 + 12H 2 = 4PH 3 + 6H 2 0. (35) 

In CH 4 -dominated objects, the PH 3 =P 4 6 equal abundance boundary is approximated by 

logP T « -11.40 + 12360/T + 2[0/H] + [P/H] 
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« -11.40 + 12360/T + 3[M/H], (36) 

assuming [O/H] ~ [P/H] (see §2). The net reaction of P4O6 to P 2 occurs at lower Pt in 
CO-dominated objects: 

P 4 6 + 6H 2 = 2P 2 + 6H 2 0, (37) 
and the position of the P 2 =P 4 6 equal abundance boundary is approximated by 

logP T « 45.17 - 51048 /T - 6[0/H] - [P/H] 

« 45.17- 51048 /T-7[M/H]. (38) 

assuming [O/H] w [P/H] (see §2). The relative importance of reactions (35) and (37) depends 
on the position of an object's P-T profile relative to the PH 3 -P 4 06-P 2 triple point at T ~ 
1101 K and Pt ~ 10 -120 bar, where all three gases have equal abundances [A(PH 3 ) = 
A(P 4 6 ) = A(P 2 ) a* ySP]. The positions of all the triple points shown in Figure 4 are listed 
in Table 2. 

In L dwarf atmospheres, PH3 and P 2 have similar abundances and are converted into 
each other by the net thermochemical reaction 

2PH 3 = P 2 + 3H 2 . (39) 

The PH 3 =P 2 equal abundance boundary in a solar metallicity gas is approximated by 

logP T « -1.56 x 10 7 /T 2 + 2.13 x 10 4 /T - 7.68, (40) 

from 1101 to 1330 K, and shifts to higher Pt and lower T with increasing metallicity. At 
deeper atmospheric levels, PH 2 becomes the most abundant phosphorus gas. In T dwarfs 
and cool L dwarfs PH 2 replaces PH 3 via 

PH 3 = PH 2 + 0.5H 2 . (41) 

The PH 3 =PH 2 boundary is given by 

log P T « 7.71 - 10888/T, (42) 

and is independent of metallicity because the m dependence for each P-bearing gas in reaction 
(41) cancels and the H 2 abundance is essentially constant over small metallicity variations. 

In Pegasi planets and hot L dwarfs, PH 2 replaces P 2 via 

P 2 + 2H 2 = 2PH 2 . (43) 
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The position of the P 2 =PH 2 boundary in a solar metallicity gas is approximated by 

logP T « -1.68 x 10 7 /T 2 + 2.62 x 10 4 /T - 10.68. (44) 

from 1330 to 1780 K, and shifts to higher P T and T with increasing metallicity. The PH 3 =P 2 , 
PH 3 =PH 2 , and P 2 =PH 2 boundaries intersect at T ~ 1330 K and P T ~ lO" 048 bar, where 
all three gases have equal abundances [A(PH 3 ) = A(PH 2 ) = A(P 2 ) \T,P]. 

At even higher temperatures, PH 2 thermally dissociates via 

PH 2 = P + H 2 , (45) 

and monatomic P becomes the major P-bearing gas. The position of the PH 2 =P boundary 
given by 

log P T « 5.00 - 1 1233/T, (46) 

is independent of metallicity because the to dependence for each P-bearing gas cancels out in 
reaction (45). As for sulfur, thermal ionization of P only becomes important at temperatures 
higher than those shown in Figure 4, e.g., P+/P ~ 1 at 6088 K (10~ 2 bars), 4836 K (10~ 4 
bars), and 4054 K (10~ 6 bars). 



4-1.1. Tetraphosphorus Hexaoxide, P^Oq gas 

The equilibrium abundances of P40 6 as a function of T and P T for solar metallicity are 
illustrated in Figure 5a. At high temperatures (> 1000 K), the P/iC) 6 abundances depend on 
which P-bearing gas is dominant. For example, within the PH 3 field, the P4O6 abundance is 
controlled by reaction (35) and is proportional to P T 3 and m 10 , assuming [O/H] ~ [P/H] (see 
§2). Inside the P 2 field, the P4O6 abundance is controlled by reaction (37) and is proportional 
to P T and m 8 , assuming [O/H] w [P/H] (see §2). In both cases, the P4O6 abundance rapidly 
decreases with temperature and is strongly dependent on metallicity. 

At low temperatures (< 1000 K), P4O6 is the dominant P-bearing gas at equilibrium 
(Figure 4), and its abundance is proportional to to, assuming [O/H] [P/H] (see §2). The 
maximum P4O6 abundance is Ap 4 o 6 ~ \X^-p because each molecule of P4O6 contains four 
phosphorus atoms. However, completely insignificant amounts of P4O6 are expected in the 
upper atmospheres of giant planets and T dwarfs because the PH 3 to P40 6 conversion is 
kinetically inhibited inside the PH 3 stability field where essentially no P40 6 exists (see §4.3). 
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4-1.2. Phosphine, PH 3 gas 

Chemical equilibrium abundances of PH3 are shown in Figure 5b as a function of T and 
Pt- At high temperatures (T > 1000 K), PH3 abundances generally increase with Pt but 
decrease with T when other P-bearing gases are dominant. For example, within the P2 field 
the PH 3 abundance is governed by reaction (39) and is approximated by 

logXpH 3 « -9.03 + 3737/T + logP T + 0.5 [P/H], (47) 

proportional to P T and m 5 . At higher total pressures, PH 3 is the dominant P-bearing gas 
(see Figure 4), and its abundance is proportional to m. The maximum PH 3 abundance is 
XpH 3 ~ A SP . Phosphine therefore approximately represents the total phosphorus inventory 
in the deep atmospheres of giant planets and T dwarfs. Furthermore, rapid vertical mixing 
and phosphine quenching from deep atmospheric regions where PH3 dominates is expected 
to give disequilibrium abundances of Aph 3 ~ A^p in the observable atmospheres of giant 
planets and T dwarfs (see §4.3). 

4-1.3. Diatomic Phosphorus, P2 gas 

Mole fraction contours for P2 are given in Figure 5c as a function of T and Pt- In cooler 
objects (e.g., giant planets and T dwarfs), P2 is governed by reaction (39), and the mole 
fraction abundance of P 2 within the PH 3 field is 

logXp 2 « -1.01 - 7474/T - 21ogP T + 2[P/H], (48) 

proportional to P^ 2 and m 2 . At lower total pressures (< 10~ ' 48 bar), P2 is the dominant 
P-bearing gas (see Figure 4), and its abundance is proportional to m. The maximum P 2 
abundance is Xp 2 m |A 2P because each molecule of P 2 contains two phosphorus atoms. 
With increasing temperatures, P 2 is replaced by PH 2 via reaction (43) in the atmospheres 
of L dwarfs and Pegasi planets. 

4.1.4. Phosphino, PH 2 gas 

Figure 5d illustrates the abundances of the PH 2 radical, which is the most abundant 
P-bearing gas in the deep atmospheres of brown dwarfs and Pegasi planets. At lower tem- 
peratures, the PH 2 abundance within the PH 3 field is governed by reaction (41) and is given 
by 

logXpH 2 « -2.40 - 5444/T - 0.5 log P T + [P/H], (49) 
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proportional to P T ' 5 and m. When P 2 is the dominant phosphorus gas, the PH 2 abundance 
is governed by reaction (43) and is given by 

logXpH 2 » -5.19 - 1707/T + 0.5 logP T + 0.5[P/H], (50) 

proportional to P? 5 and m - 5 . Equations (49) and (50) show that the PH 2 abundances 
decrease with P? when PH 3 dominates, and increase with P? when P 2 dominates. This 
explains the curvature of the PH 2 mole fraction contours in Figure 5d. 

4.2. Phosphorus Condensation Chemistry 

At low temperatures (Figure 4), equilibrium condensation of NH 4 H 2 P0 4 is expected to 
occur via the net thermochemical reaction 

10H 2 O + 4NH 3 + P 4 6 = 4NH 4 H 2 P0 4 + 4H 2 (51) 

The condensation temperature of NH 4 H 2 P0 4 as a function of P T and metallicity is approx- 
imated by 

10 4 /T cond (NH 4 H 2 PO 4 ) « 29.99 - 0.20 (11 logP T + 10[O/H] + 4[N/H] + 1[P/H]) 
10 4 /7 1 cond(NH 4 H 2 PO 4 ) « 29.99 - 0.20 (11 log P T + 15[M/H]) , (52) 

assuming [O/H] = [N/H] = [P/H] (see §2). If equilibrium were attained (which is not the 
case for the giant planets in our Solar System), NH 4 H 2 P0 4 condensation would efficiently 
remove phosphorus from the upper atmospheres of giant planets and the coolest T dwarfs 
because oxygen and nitrogen abundances are typically much larger than that of phosphorus. 
We found that neither P 4 Oio nor elemental P condense under equilibrium conditions over 
this P-T range from a solar system composition gas. 

4.3. Phosphorus Chemistry in Substellar Objects 

Figures 6a-6d illustrate equilibrium phosphorus gas chemistry along model atmosphere 
profiles Jupiter, Gliese 229B, HD209458b, and the L dwarf as done for sulfur in Figure 3. As 
before, the shape of the P-T profile of HD209458b (Figure 4) is responsible for the curvature 
in the phosphorus abundances near 1800 K. 

Using the effective temperatures for reference, the major P-bearing gases expected from 
thermochemical equilibrium in the observable atmospheres are: P 4 C>6 in giant planets and 
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cool T dwarfs (T eff ~ 100-1000 K), PH 3 in hot T dwarfs and cool L dwarfs (T eff ~ 1000-1400 
K), PH 2 in hot L dwarfs (T cff ~ 1400 - 2000 K), P 2 in Pegasi planets such as HD209458b 
(T e gf ~ 1350 K), and monatomic P in cool M dwarfs (e.g., T e g- ~ 2200 K). The dominant P- 
bearing gas changes with increasing temperature, and the major phosphorus gases at depth 
are: PH2 in Pegasi planets, PH2 and P in L dwarfs, PH3 and PH2 in T dwarfs, and PH3 in 
giant planets (Figure 4). 

However, thermochemical equilibrium does not apply to phosphorus chemistry in the 
upper atmospheres of giant planets and T dwarfs because the conversion from PH 3 to P4O6 is 
kinetically inhibited (Prinn & Owen 1976; Barshay & Lewis 1978; Prinn et al. 1984; FL94), 
and timescales for the PH3 to P4O6 conversion (t c hem) are larger than typical convective 
mixing timescales (t mix ) in these atmospheres (Fegley & Prinn 1985; FL94; Fegley & Lodders 
1996). The chemical lifetime of PH 3 as a function of T and P T for solar metallicity is given in 
Figure 7, using PH 3 destruction kinetics described by Visscher & Fegley (2005). The quench 
temperature (where t c hcm = ^mix) for PH 3 destruction is expected to be inside the PH 3 field 
in the deep atmospheres of giant planets and T dwarfs. As a result, virtually no P4O6 forms 
and PH 3 is mixed upward into their observable atmospheres at abundances approximately 
representative of the total atmospheric phosphorus inventory (Aph 3 ~ A^p). We therefore 
expect PH 3 (instead of P4O6) to be the major P-bearing gas in giant planets and T dwarfs. 

Indeed, spectroscopic observations of PH 3 in the upper atmospheres of Jupiter and 
Saturn indicate tropospheric phosphine abundances of 1.1 ppm by volume on Jupiter (e.g., 
Lodders 2004b, and references therein) and 4.5 ppm by volume on Saturn (e.g., Visscher & 
Fegley 2005, and references therein). These measured PH 3 abundances are ~ 30 orders of 
magnitude higher than expected from thermodynamic equilibrium and imply rapid vertical 
mixing from deeper atmospheric levels where phosphine is more abundant. We expect PH 3 
abundances of ~ 0.6 ppm (for solar metallicity) in the upper atmospheres of Jupiter-type 
giant planets and cool T dwarfs such as Gliese 229B (e.g., Fegley & Lodders 1996), and 
measurements of phosphine in these objects should approximately represent their bulk at- 
mospheric phosphorus inventory. Noll & Marley (1997) discussed the observability of gases 
affected by vertical mixing and concluded that the 4.3 /im feature of PH 3 is potentially de- 
tectable in the infrared spectra of T dwarfs. Photochemistry may plausibly destroy PH 3 in 
the uppermost atmospheres of these objects (see §5), as occurs on Jupiter and Saturn. 

In hotter objects, PH 3 may no longer represent the total atmospheric phosphorus in- 
ventory when P equilibrium chemistry is established and other abundant P-bearing gases 
are present. Phosphine abundances of ~ 50 ppb (for solar metallicity) remain at the 1350 
K level on HD209458b and 1800 K in the model L dwarf. The PH 3 feature at 4.3 /im is 
potentially detectable in these objects if it can be distinguished from nearby CO bands (Noll 
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& Marley 1997; Saumon et al. 2003). 

The abundant P-bearing gases at higher atmospheric temperatures include P 2 , PH, PH 2 , 
PH 3 , and HCP (Figures 6a-6d). The phosphino radical (PH 2 ) is the dominant P-bearing 
gas in the observable atmosphere of the model L dwarf and in the deep atmospheres of T 
dwarfs and Pegasi planets. However, PH 2 only represents about half of the total atmospheric 
phosphorus inventory because other P-bearing gases remain relatively abundant within the 
PH 2 field. Equilibrium PH 2 abundances of ~ 1 ppb, ~ 50 ppb, and ~ 300 ppb are expected 
(for solar metallicity) in the observable atmospheres of T dwarfs, Pegasi planets, and L 
dwarfs, respectively. 

As temperatures decrease, PH 2 is replaced by PH 3 (in Gliese 229B, Figure 6b) and/or P 2 
(in HD209458b and the L dwarf, Figures 6c-6d) as the major P-bearing gas, depending on the 
position of the atmospheric P-T profile relative to the PH 3 -PH 2 -P 2 triple point (see Figure 
4 and §4.1). Although P 2 is expected to be the dominant P-bearing gas in the observable 
atmospheres of Pegasi planets (Xp 2 ~ 0.2 ppm at 1350 K) and the upper atmospheres of L 
dwarfs, the P 2 bands expected near 0.2 /zm (Carroll & Mitchell 1975) may be very difficult 
to detect in these objects. At temperatures below ~ 1000 K, P 2 is converted into P40 6 , 
which is plausibly destroyed by photochemical processes in the uppermost atmospheres of L 
dwarfs and Pegasi planets (see §5). 

The abundances of PH and HCP are negligible in the observable atmospheres of giant 
planets and T dwarfs. In warmer objects, these gases have equilibrium abundances (at solar 
metallicity) of ~ 2 ppb (PH) and ~ 20 ppb (HCP) on HD209458b (at 1350 K), and ~ 50 
ppb (PH) and ~ 15 ppb (HCP) in the L dwarf model atmosphere (at 1800 K). 

Under equilibrium conditions, NH 4 H 2 P0 4 cloud formation would be expected near the 
400 K level in the upper atmosphere of Jupiter (cf. Barshay & Lewis 1978; FL94), but is 
precluded by quenching of PH 3 (see above). Condensation of NH4H 2 P04 is not expected 
in the warmer atmospheres of Gliese 229B, HD209458b, and the L dwarf because their P-T 
profiles do not intersect the NH 4 H 2 P04 condensation curve (Figure 4). 

5. Thermochemistry and Photochemistry 

So far, our treatment of sulfur and phosphorus chemistry assumes the establishment 
of thermochemical equilibrium. The exception is phosphine, for which the kinetics of its 
destruction mechanism and convective vertical mixing must be considered for determining 
PH 3 abundances (§4.3). However, we must also consider photochemical reactions in the 
uppermost atmospheres of planets and brown dwarfs that are companions to stars. If ul- 
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traviolet flux causes photochemical reactions, the thermochemical reactions will be driven 
out of equilibrium. Photochemistry may also play a role in isolated brown dwarfs if their 
coronae are warm enough to produce an ultraviolet flux (Yelle 1999). 

For example, in Jupiter's troposphere the major C-, N-, S- and O-bearing gases ex- 
pected from thermochemical equilibrium are CH 4 , NH 3 , H 2 S, and H 2 0, respectively (FL94). 
However, Jupiter's stratosphere contains photochemically produced hydrocarbons such as 
C2H 6 , C2H2, and C2H4 (Gladstone et al. 1996). Hence we must examine the relative roles of 
thermochemistry and photochemistry on S and P compounds in brown dwarfs, EGPs, and 
related objects. 



5.1. Thermochemistry in Substellar Objects 

Thermochemical equilibrium is readily achieved in the hot, deep atmospheres of sub- 
stellar objects because gas phase reaction rates are generally much faster than convective 
mixing rates (i.e., t c h em <C t m i X ). To illustrate this, we consider the kinetics of the following 
elementary reactions: 

H + H 2 S -> SH + H 2 , (53) 
SH + CO -> OCS + H, (54) 
S + H 2 S -> SH + SH, (55) 
Na + HCl -> NaCl + H, (56) 
PH + PH 2 -> PH 3 + P, (57) 
H + PH 3 -> PH 2 + H 2 . (58) 

Here we use reaction (53) as an example for computing t c h em - The rate of disappearance of 
H 2 S in reaction (53) is given by 

d[H 2 S] 



dt 



-MHpaS], (59) 



where (cm 3 molecules -1 s -1 ) is the rate coefficient for reaction (53) and [H 2 S] and [H] 
are the molecular number densities (molecules cm -3 ) for H 2 S and H. The chemical lifetime 
for H 2 S is defined as 

WH 2 S) = - d[Ha g ]M = T^pj, (60) 

and is the time required for a 1/e reduction in the H 2 S abundance. Representative chemical 
lifetimes for reactions (53)-(58) at 1 bar total pressure and 1000, 1500, and 2000 K are listed 
in Table 1. Although the short chemical lifetimes in Table 1 indicate that these reactions 
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proceed rapidly, we must also consider vertical mixing rates to determine if equilibrium is 
established. The timescale for convective mixing, t mix is approximated by 

tmix ~ H 2 1 'K eddy (61) 

where K e dd y is the eddy diffusion coefficient. The pressure scale height, H, is given by 

H = RT/fig (62) 

where R is the gas constant, /x is the mean molecular weight of the gas, and g is gravity. 
For a solar system composition gas, fx xs 2.3 g mol -1 . A convective timescale of t mix ~ 10 6 
seconds is expected near the 1000 K level on Jupiter, using H ~ 140 km (logg = 3.4) and 
Keddy ~ 10 8 cm 2 s _1 (Fegley & Prinn 1988). On Gliese 229B, we estimate t m i X ~ 10 7 seconds 
at the 1000 K level, assuming H ~ 3.6 km (log g = 5.0; Saumon et al. 2000) and K e dd y ~ 10 4 
cm 2 s _1 (Griffith & Yelle 1999). Because of its low gravity (logg = 3.0; Iro et al. 2005) 
HD209458b has large mixing times of ~ 10 7 to 10 11 seconds at the 1000 K level, depending 
on the assumed K eddy value (~ 10 8 to 10 4 cm 2 s _1 ). These timescales for convective mixing 
are typically many orders of magnitude greater than the chemical lifetimes listed in Table 
3 (i.e., t c hem *C t mix ). We therefore expect the rapid establishment of t her mo chemical 
equilibrium in the deep atmospheres of substellar objects. In instances where the rate of 
vertical mixing equals or exceeds the rate of chemical reactions (t c h em > t mix ), disequilibrium 
species (e.g., PH3, see §4.3) are mixed into the upper atmospheres of substellar objects, as 
occurs on Jupiter and Saturn (e.g., CO, HCN, PH 3 , GeH 4 , AsH 3 ; FL94) and Gliese 229B 
(e.g. CO, PH 3 , N 2 ; Fegley & Lodders 1996; Griffith & Yelle 1999; Saumon et al. 2000). 



5.2. Photochemistry on Pegasi Planets 

In order to examine the relative importance of photochemical vs. thermochemical pro- 
cesses in substellar objects, we consider an extreme example: the strongly irradiated upper 
atmosphere of HD209458b. This Pegasi planet orbits a solar-type star at 0.05 AU and 
receives a stellar UV flux ~ 10,000 times that for Jupiter (e.g., Liang et al. 2003). To de- 
termine the depth in the atmosphere where thermochemical processes become important, 
we compared the thermochemical equilibrium abundances of the reactive species OH, O, 
and H in our model with photochemical equilibrium abundances from Liang et al. (2003) 
for HD209458b. Plots comparing photochemical (dotted lines) and thermochemical (solid 
lines) abundances for these reactive species are shown in Figure 8, using the P-T profile for 
HD209458b from Figure 1 of Liang et al. (2003). 

At very low pressures, in the uppermost atmosphere of HD209458b, photochemistry 
produces excess disequilibrium abundances of OH, O, and H, primarily via H 2 photoly- 
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sis (Liang et al. 2003). Thermochemistry becomes increasingly important with depth, as 
pressures and temperatures increase and the available stellar flux decreases. Eventually, the 
predicted thermochemical abundances surpass the photochemical abundances, indicating the 
establishment of equilibrium chemistry. In other words, the photochemical reactions are no 
longer driving the thermochemical reactions out of equilibrium. As shown in Figure 8, this 
crossover occurs near the ~ 1 to 10 mbar levels (T ~ 980 to 1140 K) for the reactive species 
OH, O, and H. Thermochemical processes are therefore expected to be important at pres- 
sures of ~ 10 mbar and higher on HD209458b, consistent with results showing that most of 
the incoming stellar flux is absorbed by the 1 bar level (T ~ 1640 K) in its atmosphere (Iro 
et al. 2005). Furthermore, Fortney et al. (2005, 2006) point out that the infrared spectra of 
Pegasi planets appear to be most sensitive to opacity near 10 to 100 mbar (T ~ 1140 to 1400 
K). This suggests that thermochemical models will provide a useful basis for interpreting 
and guiding observations of Pegasi planet atmospheres. 



6. Summary 

Thermochemistry governs the chemical behavior of sulfur and phosphorus species in the 
hot, deep atmospheres of substellar objects. Hydrogen sulfide is the most abundant sulfur- 
bearing gas throughout substellar atmospheres, and observations of H 2 S in these objects 
should provide a good estimate of their atmospheric sulfur content. The condensation of 
metal sulfide clouds slightly lowers the atmospheric H 2 S abundance to ~ 91% of the total 
S inventory. These clouds also affect opacity by introducing condensed particles and by 
removing absorbing gases from the observable atmosphere. Mercapto (SH) and SiS are the 
next most abundant sulfur gases after H 2 S, and the relative importance of SiS increases with 
increasing effective temperature. Furthermore, the maximum SiS abundance occurs near the 
Mg-silicate cloud base and is a potential tracer of weather in Pegasi planets and L dwarfs. 

Phosphorus speciation is considerably more complex than that for C, N, O, or S, and 
several phosphorus gases become relatively abundant at high temperatures. Disequilibrium 
abundances of PH 3 are expected in the upper atmospheres of giant planets and T dwarfs 
via rapid vertical mixing from deeper levels. In addition, phosphine abundances in giant 
planets and T dwarfs are expected to approximately represent their atmospheric phosphorus 
inventories. In Pegasi planets and L dwarfs, P 2 is the dominant phosphorus gas until replaced 
by P4O6 at low temperatures or PH 2 at high temperatures. Phosphino (PH 2 ) is the most 
abundant phosphorus gas in the deep atmospheres of brown dwarfs and Pegasi planets. 

This work was supported by the NASA Planetary Atmospheres Program. Work by KL 
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A. Additional Compounds 

Here we list additional sulfur and phosphorus compounds investigated but not discussed 
in detail above. This list includes unstable condensates and minor gases predicted to have 
very low abundances in a solar system composition gas at temperatures and total pressures 
relevant for substellar atmospheres. Depending on the prevailing sulfur, carbon, and phos- 
phorus chemistry, the abundances of key gases are Pt independent: H2S2 within the H2S 
field; OCS, CS, and CS 2 within the CO field; CH 3 SH within the CH 4 field; HPO, PH, and 
PN within the P2 field. However, their abundances are generally too low to serve as useful 
temperature probes in brown dwarf or EGP atmospheres. Some gases (H2S 2 , CH3SH) are 
potentially detectable via deep atmospheric entry probes on Jupiter and Saturn (FL94). 

For a given metallicity, relatively constant OCS abundances are expected throughout 
the atmospheres of CO-dominated objects; however, maximum OCS abundances are only 
~ 1 ppb (for solar metallicity) in L dwarfs and Pegasi planets. The equilibrium abundances 
of P, PO, and PS generally increase with T and decrease with Pt- These gases only become 
important at lowest total pressures considered here (see Figure 4). The behavior of HCP is 
similar to that of PH 2 ; maximum HCP abundances are ~ 50 ppb (for solar metallicity) over 
the P-T range considered here. 

Sulfur 

S 3 , S 4 , S 2 0, SO, S0 2 , S0 3 , SOH, HSO, H 2 SO, H 2 S0 4 , SF, SF 2 , H 2 S 2 , OCS, CS, CS 2 , CH 3 SH, 
NS, PS, K 2 S, CaS, FeS, MgS, MnS, TiS, VS, ZrS, SiS 2 , TiS 2 , ZrS 2 , S(s,l), OCS(s,l), S0 2 (s,l) 

Phosphorus 

P 3 , P 4 , PO, P0 2 , P2O3, P 2 4 , P 2 5 , P 3 6 , P 4 7 , P 4 8 , P 4 9 , P 4 O 10 , PCI, PC1 2 , PCI3, 
PF, PF 2 , PF 3 , HPO, HCP, CP, PH, PN, PS, P(s), P 4 Oi (s,l) 
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Table 1. Equilibrium Condensation Temperatures of Sulfur-Bearing Compounds 



object MnS Na 2 S ZnS NH 4 SH 



Jupiter 1800 

Gliese 229B 1380 

HD209458b 1260 

L dwarf 1270 
M dwarf 



1370 970 220 
1000 780 
~ 800 a 



Note. - - Condensation temperatures com- 
puted using [M/H] « +0.5 for Jupiter, [M/H] 
« -0.3 for Gliese 229B, and [M/H] « for 
HD209458b and the model L dwarf. 

a night-side condensation 



Table 2. Phosphorus Equal Abundance "Triple Points" in a Solar Metallicity Gas 



triple point a 


T, K 


log(P T ,bar) 


PH 3 -P 4 6 -P 2 


1101 


-1.20 


PH 3 -PH 2 -P 2 


1330 


-0.48 


P 2 -PO-PS 


1389 


-3.93 


P 2 -PH 2 -PO 


1780 


-1.26 


PH 2 -P-PO 


1811 


-1.20 



a e.g., where Xa = X& = Xq for 
phosphorus-bearing gases A, B, C. 
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Table 3. Representative Chemical Lifetimes at 1 bar pressure 



chemical lifetime, s 



reaction 


tchem 


1000 K 


1500 K 


2000 K 


kinetic data 


53 


tchem (H2S) 


1Q 0.4 


1Q -3.8 


10 -5.9 


Yoshimura et al. (1992) 


54 


tchem (SH) 


io 2 - 7 


1Q 0.2 


1Q -0.2 


Kurbanov & Mamedov (1995) 


55 


^c/iem(S) 


10 L7 


1Q -2.3 


1Q -4.3 


Woiki & Roth (1994) 


56 


tchem (Na) 


io 1 - 4 


10 -1.6 


10 -1.8 


Husain & Marshall (1986) 


57 


tchem (PH) 


io 4 - 2 


io 1 - 2 


IO 10 


Twarowski (1995) 


58 


^chem(PH 3 ) 


10 0.4 


1Q -3.5 


1Q -5.4 


Arthur & Cooper (1997) 



Note. — Convective mixing timescales (in seconds) are t conv ~ 10 6 ' 3 (1000 K), 10 6 ' 6 
(1500 K), and IO 69 (2000 K) for Jupiter, 10 71 (1000 K), 10 7 - 5 (1500 K), and 10 7 - 7 
(2000 K) for Gliese 229B, and 10 7 - 1 to 10 11 ' 1 (1000 K), 10 7 - 5 to IO 115 (1500 K), and 
10 7 - 7 to 10 11 - 7 (2000 K) for HD209458b. See §5.1 for details. 
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log (P T , bar) 



Fig. 1. — Overview of sulfur chemistry as a function of temperature and pressure in a solar 
system composition gas. The solid lines indicate where major S-bearing gases have equal 
abundances. The triangle denotes the position of the triple point where H 2 S, SH, and S have 
equal abundances. Also shown are condensation curves for S-bearing compounds (dotted 
lines) and the H 2 =H and CH 4 =CO equal abundance lines (dash-dot lines). Atmospheric 
profiles (dashed lines) are given for Jupiter (T eff = 124 K), Gliese 229B (T cff = 960 K), 
HD209458b (T cfr = 1350 K), an L dwarf (T cS = 1800 K), and an M dwarf (T cff = 2200 K). 
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Fig. 2. — Mole fraction contours (on a logarithmic scale) for (a) hydrogen sulfide (H 2 S), (b) 
mercapto radical (SH), (c) monatomic sulfur (S), (d) silicon sulfide (SiS) as a function of 
pressure and temperature in a solar metallicity gas. The dashed lines labeled Fo and En in 
the SiS plot indicate the equilibrium condensation temperatures of forsterite (Mg 2 Si0 4 ) and 
enstatite (MgSi0 3 ), respectively. 
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Fig. 3. — Sulfur chemistry along representative atmospheric pressure-temperature profiles 
from 900 to 2500 K for (a) Jupiter ([M/H] w +0.5), (b) Gliese 229B ([M/H] w -0.3), (c) 
HD209458b ([M/H] w 0), and (d) an L dwarf ([M/H] w 0; T cfr = 1800 K), assuming uniform 
elemental enrichments and depletions. Also shown are the calculated condensation tem- 
peratures of forsterite (Fo; Mg 2 Si0 4 ), enstatite (En; MgSi0 3 ), MnS, Na 2 S, and ZnS. With 
the exception of forsterite (not shown) and enstatite on Jupiter, all condensates shown here 
condense as solids. 
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NH 4 H 2 P0 4 (s) - 




log {P v bar) 



Fig. 4. — Overview of phosphorus chemistry as a function of temperature and pressure in 
a solar system composition gas. The solid lines indicate where major P-bearing gases have 
equal abundances. Also shown is the condensation curve for NH 4 H 2 P0 4 (dotted line) and the 
position of the H 2 =H equal abundance boundary (dash-dot line). Atmospheric profiles for 
representative substellar objects are shown for reference (dashed lines). See text for details. 
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Fig. 5. — Mole fraction contours (on a logarithmic scale) for (a) tetraphosphorus hexaoxide 
(P4O6), (b) phosphine (PH 3 ), (c) diatomic phosphorus (P2), and (d) phosphino (PH 2 ) as a 
function of pressure and temperature in a solar metallicity gas. Gas abundances at higher 
or lower metallicities are given by the abundance expressions and/or the dependence on the 
metallicity factor m. 
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Fig. 6. — Phosphorus equilibrium chemistry along model atmosphere P-T profiles from 800 
to 2500 K for (a) Jupiter ([M/H] « +0.5; T cff = 124 K), (b) Gliese 229B ([M/H] « -0.3; T eff 
= 960 K), (c) HD209458b ([M/H] « 0; T eff = 1350 K), and (d) an L dwarf ([M/H] « 0; T eff 
= 1800 K), assuming uniform elemental enrichments and depletions. Note that the PH3 to 
P4O6 conversion is kinetically inhibited in giant planets and T dwarfs, and that Jupiter's 
observable atmosphere lies outside the temperature range shown here. 



-37- 
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Fig. 7. — Logarithmic time scale (seconds) for the chemical conversion of PH3 to P4O6 in 
a solar metallicity gas. This conversion is kinetically inhibited in the upper atmospheres of 
giant planets and T dwarfs. 
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Fig. 8. — Comparison of thermochemical (solid lines) and photochemical (dotted lines; Liang 
et al. 2003) abundances for the reactive species OH, O, and H in the upper atmosphere of 
HD209458b. See text for details. 



